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Abstract. Gaussian graphical models are parametric statistical models for jointly nor- 
mal random variables whose dependence structure is determined by a graph. In previous 
work, we introduced trek separation, which gives a necessary and sufficient condition in 
terms of the graph for when a subdeterminant is zero for all covariance matrices that 
belong to the Gaussian graphical model. Here we extend this result to give explicit 
cancellation-free formulas for the expansions of nonzero subdctcrminants. 



1. Introduction 

Gaussian graphical models are parametric statistical models for jointly normal random 
variables whose dependence structure is determined by a graph. In this work, we con- 
sider Gaussian graphical models on mixed graphs which have both directed edges and 
bidirected edges. In the causal modeling framework for graphical models [H], the directed 
edges represent direct causal effects of one variable on another, while the bidirected edges 
represent the effects of unobserved confounders that lead to correlations between the er- 
ror terms of the variables. In other contexts, the Gaussian graphical models we study 
are known as linear structural equation models, and the directed edges describe linear 
relationships between variables with correlated errors. 

A fundamental problem in the study of graphical models is to characterize the distri- 
butions that can arise. The models are typically presented using parametric descriptions, 
and we wish to give conditions on the distributions which could appear for some choice 
of parameters. In the case of graphical models which only have directed edges without 
directed cycles (so-called directed acyclic graphs or DAGs), it is well-known that a prob- 
ability distribution belongs to the model if and only if the distribution satisfies all the 
conditional independence constraints implied by the graph [G*, Thm 3.27]. 

For Gaussian graphical models, the model consists of all covariance matrices S G PDm 
which arise for some choice of parameters. Each conditional independence statement 
X^iLX^lXc, translates to a rank condition on a submatrix of the covariance matrix E, 
namely that the matrix ^auc,buc has rank < i^C. (Note that for sets I and J, S/^j 
denotes the submatrix of S with row index set I and column index set J.) In other 
words, all (#C + 1) x (#C -|- 1) subdeterminants of T,auc,buc are zero. 

In previous work [9j, we generalized the rank dropping condition from conditional in- 
dependence to arbitrary minors, via trek separation. This gives necessary and sufficient 
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conditions in terms of the underlying graph G, and sets A and B, such that det T,a,b = 0. 
These vanishing determinantal constraints are among the most natural to study. In the 
present paper, we extend the result of [U] to give explicit cancellation-free expansions for 
the no n- vanishing determinants det Tia,b- This has potential applications to the algebraic 
study of graphical models including to the identifiability problem for Gaussian graphical 
models [H |3] . Besides this, our formulas for det T,a,b involve elegant combinatorics in the 
spirit of classical enumerative combinatorics results on determinants of matrices associ- 
ated to graphs [U [7]. For example, our results give the expansion of det ^ in terms 
of nonintersecting trek flows, and we show that after collecting terms, every monomial 
in the expansion appears with a coefficient of the form ±2'^ for some c, determined by 
combinatorial properties of the flow. 

This paper is organized as follows. In the next section, we give an introduction to 
Gaussian graphical models and the special "paths" in mixed graphs we need to study, 
called treks. In Section [3] we describe our main results on the expansions of the deter- 
minants of submatrices of S arising from a Gaussian graphical model. These results are 
based on a generalization of the Lindstrom-Gessel-Viennot Lemma for graphs with cycles, 
which is described in Section |4j Finally, Section |5] contains the proofs of the most general 
statements from Section |3] using the methods of Section |4j 

2. Gaussian Graphical Models 

Let G = {V,B,D) be a mixed graph with vertex set V, bidirected edge set B, and 
directed edge set D. Thus (V, D) is a directed graph without loops or multiple edges 
called the directed part of G, and {V, B) is an undirected graph without loops or multiple 
edges called the bidirected part of G. The qualiflcation bidirected stems from these edges' 
statistical interpretation in graphical models. A bidirected edge between vertices i and j 
is denoted ^ O j and a directed edge from i to j is denoted i ^ j. Let PDy denote the 
set of X j^V real symmetric positive deflnite matrices. Let PD{B) denote the set of 
positive deflnite matrices with off-diagonal nonzero entries only in positions corresponding 
to bidirected edges of G, i.e. 

PD{B) = {Q = {uij) G PDv : = if i ^ j and z ^ j ^ B}. 

Let e G be a jointly normal random vector e ~ A/'(0, VL) with Vt G PD{B). For each 
z — )■ j G -D, let Xij G M be a real parameter. Deflne a new random vector X G MY by 

(1) X,= Yl >^^jX, + e,. 

Let A be the j^V x jj=-V matrix with Ajj in the ij position if i — )■ j G -D, and zero otherwise. 
Equation ([T]) can be rearranged as e = (/ — A)^X, where / is a jj^V x jj^V identity matrix. 
From this, we deduce that X ~ A/'(0, S) where 

S = (/- A)-^fi(/- A)-i 

provided that the matrix / — A is invertible. 

Given A,B C V, let T,a,b = {o'a,b)aeA,beB be the submatrix of S with row index 
set A and column index set B. Suppose #A = ^B. In previous work [_9J, we gave a 
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combinatorial condition for when det Ti^^b = 0. In this paper, we explain how to expand 
det TiA,B explicitly in monomials without cancellation. Henceforth, rather than thinking 
about the Ajj and Uij as parameters, we choose to think about them as indeterminates (or 
polynomial variables). Hence, we are interested in expanding det S^^^ as a polynomial or 
a rational function, depending on the context. 

One of the key combinatorial definitions we need is the definition of a trek. 

Definition 2.1. A trek in G from i^Vtoj&Visa. pair [Pl,Pji), where Pl is a 
directed path from some vertex s to i and Pr is a directed path from some vertex t to 
j, satisfying the further requirement that ii s ^ t, then there must be a bidirected edge 
s ^ t E B. The vertices i and j are called the initial vertex and the final vertex of the 
trek, respectively. 

In this definition, a path is an ordered sequence of edges in D in which the head of 
each edge equals the tail of the next edge. A path may have no edges, in which case, we 
also specify where the path starts, which is also where it ends; we call such a path an 
empty path. Later in the paper we will work mostly with self-avoiding paths, in which 
each vertex is the head of at most one edge of the path and the tail of at most one edge 
of the path; such a path is uniquely determined by its edge set (and initial vertex when 
empty). For the moment, however, Pl and Pr need not be self-avoiding. Note also that 
either (or both) of the paths Pl and Pr are allowed to be an empty path, in which case 
s = i or t = j, respectively (or both). Two treks T = {Pl, Pr) and T' = {P[, P^) are the 
same if both Pl = P'l and Pr = P'p^. In particular, two treks might contains exactly the 
same set of edges (even with the same multiplicities) but be different treks. 

If all directed edges in G are pointing down the page — which they can be arranged to 
do if the directed part of G happens to be acyclic — then one can think of a trek between 
i and j as a path of edges that starts at i, goes up for a time to s, possibly traverses a 
single bidirected edge to a vertex t, and then turns around and goes down for a time to 
j; see Figure [T] below. We stress, however, that our results below are not restricted to the 
case where the directed part of G is acyclic. 

To any trek T = [Pl,Pr), we associate the trek monomial rriT = X^Ust^^, where 
= rifc^iePi similarly, A^ = Ylk^iePn Here s,t are the initial vertices of 

Pl, Pr, and a variable Xki appears in A^ and A^ with its respective multiplicities in the 
paths Pl and Pr. The reason for introducing treks and trek monomials is the following 
result: 



Proposition 2.2. Let G = {V,B,D) be a graph, and S = (/ - A)-'^Q{I - A)-^ Then 
for all i,j G V , 




T 



where the sum runs over all treks T in G from i to j. 



Proposition 2.2 can be proven use the expansion 



(/-A)-i = J + A + A' + -- - , 
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Figure 1. Tailswapping at {k,l) yields a trek with the same trek mono- 
mial. The two treks on the right are the same. The number of treks from i 
to j with this trek system monomial equals 2**^-^^"'^^"^^ = 2^~^ = 4. 

and the combinatorial interpretation of the kl entry of A*" as the sum of A"^ over all directed 
paths P of length r from k io I ui the graph (V, D). See P]. 

As stated previously, there might be many treks in T{i,j) that use exactly the same set 
of edges, so the expression for aij might have repeated terms. As a prelude to our main 
results, we discuss how to simplify this expression in the case where (V, D) is acyclic. Let 
T = {Pl, Pr) be a trek with Pl = [s = iq ^ ii ^ . . . ^ im = i) and Pr = (t = jo — j- 
J 1 — 7- . . . — 7- j„ = j). Since iV,D) is acyclic, the ik are mutually distinct, and so are the 
jl. But it may happen that equals ji, and then T' := {P'j^,P'^) with P'j^ : t = j^ ^ 
... ^ ji-^ ik+i ... ^ im = i and P'j^ : s = io ^ . . . ^ ik ^ ji+i . . . -> j„ = j is a 
trek with m^' = mT- We call this procedure tailswapping at see Figure [l] It turns 

out that any trek T' with mx = mrpi can be obtained from T by repeated tailswapping. 
We call such treks equivalent to T, and we write T T' . \{ k = I = Q (so that s equals 
t), then then clearly T = T'. Similarly, if /c, / > 1 and ik-i = ji-i, then tailswapping at 
(A; — 1, / — 1) gives the same result as tailswapping at {k, I). Denote by i{T) the number 
of pairs {k, I) ^ (0, 0) with i^ = ji and by e(T) the number of edges that Pl and Pr have 
in common. The above discussion, together with the observation that tailswapping at 
distinct pairs of indices commutes then yields the following proposition. 



Proposition 2.3. Let G = {V, B, D) be a mixed graph such that [V, D) is directed acyclic, 
and let S = (/ - A)-^n{I - A)-\ Then for all i,j e V, 




[T] 



where the sum runs over equivalence classes of treks. 



We will greatly generalize Proposition 2^ in Section [3j In particular, the fact that all 
coefficients in det T,a^b are powers of 2 persists when (V, D) is a directed acyclic graph. 
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Now we return to the general case, where {V, D) may contain cycles. Since each entry cjjj 
is the sum of trek monomials from i to j, the determinant det Syi,s will be a (signed) sum 
of products of trek monomials. Indeed, expanding the determinant using the standard 
Leibniz expansion yields the following: 

Proposition 2.4. Let A,B C V, with — — k. Let A — {ai,...,afe} and 
B = Then 

det Ea_b = 2Z X] ^^Sn (tt) ■ uit^ ■ ■ ■ uit^ 

There are two issues with this expansion. The first is that distinct tuples (tt, Ti, . . . , T^) 
can lead to the same monomial rriTi ■ ■ ■ rriT^. and thus cancellation can occur. The second 
is that when the directed part of G contains directed cycles, this expansion is a formal 
power series rather than a polynomial (of course, this already happens when A and B 
are singletons). In the general case our aim is to write this formal power series as a 
rational function in which the denominator and the numerator are both in a suitable 
sense cancellation-free. 

An important simplification that we can make when describing the graphs, matrices, 
and determinants is to eliminate the bidircctcd edges under consideration. Associated to 
the graph G = {V. B, D), we introduce the bidircctcd subdivision G, which has only di- 
rected edges. We start will all vertices and directed edges in G. Then, for each bidirected 
edge i •<->■ J in G, we introduce a new vertex (i, j), and two directed edges (i, j) — >■ i and 
— > j. The resulting graph is the bidirected subdivision of G. Combinatorial expres- 
sions for the entries dij associated to the graph G can be obtained from the corresponding 
expressions for G as described in the following proposition. 

Proposition 2.5. Let G = {V,B,D) be a graph and let G be its bidirected subdivision. 
Write S = (/ — A)^^Q{I — A)~^ for the covariance matrix for the graph G and S* = 
(/ — A*)~-^f2*(/ — A*)"-*^ for the covariance matrix for the bidirected subdivision G. Let 
f e R[(t] be a polynomial. Then f{a{X,oj)) — if and only if f{a*{X* ,oj*)) — 0. 

Furthermore, the expansion of f{a{X,uj)) can be recovered from the expansion off{a*{X*, 
by the following procedure: 

(1) Remove any monomial that contains a (^(jj) or 

(2) Set all remaining X*^^-^. and X*^^--^- variables to 1. 

(3) Change ul-j^^^^^ .^ touij. 

(4) Change all remaining X*j and uj*^ to Xij and cuu, respectively. 

Proof. Consider the M-algebra homomorpliism (f) : ]R[[A, a;]] — > ]R[[A*, cu*]] defined by 

j-.i^j&B 

We claim that a*j = (l){(Jij). 

First, we prove the claim. Let A* and Q,* be the matrices for the bidirected subdivision 
G. We reahze both matrices as block matrices with the first set of row/column labels 
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corresponding to bidirected edges i ^ j, and hence the vertices in G. Then A* and 
Q* have the form 

^ =(o aJ ^ =(o n;) 

where fi* = diag{{u(ij)^(^ij))i^ji^B), ^2 = diag{{uJii)i(zy), A is the directed edge matrix 
from the graph G, and is a matrix whose {i,j),k entry is )^{ij),k if ^ € and is 

zero otherwise. Note that 

from which we deduce that 

(J-Aj U[I-A) = _ A)-^^^^]* {I-A)-^{Ql + EmiE){I-A)-^J- 

The entry in the bottom right-hand corner of this matrix is the expression for the sub- 
matrix Syy, so the claim is equivalent to saying that the map is that map that takes 

in (/ — Ay^Q{I — A)~^ and replaces it with F = + E'^Vt\E. But it is easy to see 
that = (piujii) and for i 7^ j, = if i j ^ 5 and 7,^ = (pioJij) if i -H- j E B. 

Now we will show that the claim proves the proposition. First of all, the map (p is 
injective. Indeed, suppose that 0(/) = and / is irreducible. Then, since a;*j only 
appears in (p^uu), then uu could not appear in /. Similarly, we can then rule out Uij 
appearing in /, and Xij appearing in /. This forces / to be the zero polynomial. Since 
is injective, we deduce the first part of the proposition, a power series f{cr{X,u)) is zero 
if and only if = f{a*{X* ,uj*)) is zero. 

For the second part, we note that the following map ip applied to f{a*{X*,uj*)) produces 
f{a{X,u)y. 

That can be seen by applying the operations on the factorization (J— A)~^(f22+-^^^i-^)(-^~ 
A)~^ to produce (J — A)~'^Q{I — A)~^. This corresponds to the four-step procedure from 
the statement of the proposition, since the expressions involving j) j) (i j) "wiU 
cancel with terms coming from ^'(w*;) = oJu — Ylij-i^j^B^ij- '-' 

Henceforth, we focus solely on the case of graphs with no bidirected edges. 



3. Results 

In this section, we assume G = (V, D) is a directed graph, possibly obtained as the bidi- 
rected subdivision of some mixed graph. Our goal is to find a cancellation-free analogue 



of the formula for detS^^ in Proposition 2.4 The exact meaning of cancellation-free 



depends on the context. If G is acyclic, then det Tja,b is a polynomial in the variables Ua 
and the variables X^j for directed edges i ^ j, and cancellation-free means that we deter- 
mine which monomials have non-zero coefficients in det ^a,b, and what those coefficients 
are. 
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If, on the other hand, G contains directed cycles, then there may be infinitely many 
treks from i to j. In this case det ^ is a formal power series representing a rational 
function (as is clear from E = (J — A)^'^Q{I — A)^^ using Cramer's rule). We will then 
write det T,a,b as a fraction ^ where R and S are polynomials with known monomials 
and known coefficients. 

The formula that we derive in the general case has 5 = 1 when specialized to the 
acyclic case, hence our proof will immediately focus on the general case. For the purpose 
of exposition, however, we first present the formula for the acyclic case, before dealing 
with the general case, which is more combinatorially complicated. 

3.1. Directed acyclic graphs. 

Definition 3.1. Let A and B be sets of k vertices. A trek system T from AtoB consists 
of k treks whose initial vertices exhaust the set A and whose final vertices exhaust the set 
B. 

Such a trek system gives rise to a bijection A B mapping a to the final vertex of 
the (unique) trek in T starting at a. Given a linear ordering oi, . . . , of A and a linear 
ordering bi, . . . ,bk of B, this bijection determines an element vr of &k (which depends 
on the linear orderings of A and B), and we define signT := signvr. We define the trek 
system monomial rriT as the the product of the trek monomials rriT ranging over T G T. 



The formula in Proposition |2.4| can then be rewritten as 
(2) det S^,B = ^ sign (T)mT 

T 

where T runs over all trek systems from A to B. Cancellation can happen if distinct mx 
and mx' are the same. To achieve a cancellation- free formula we introduce the following 
notion (see [9]). 

Definition 3.2. A trek system T from A to B has no sided intersection if the left 
parts Pl for (P^, Pr) G T are mutually vertex-disjoint and also the right parts Pr for 
{Pl, Pr) G T are mutually vertex-disjoint (but any Pl may have vertices in common with 
any P^). Otherwise, T is said to have a sided intersection. We write T{A, B) for the set 
of trek systems from A to B with no sided intersection. 

In [9] it is proved that in ^ one may restrict T to run over all trek systems without 
sided intersection. Our first new result is that these terms all have the same sign and 
thus do not cancel. 

Theorem 3.3 (Positivity for acyclic digraphs). Let G = (y,D) be an acyclic digraph and 
let A and B be subsets ofV of the same cardinality. Assume we have fixed linear orderings 
of A and B. If two trek systems T and T' from A to B without sided intersection satisfy 
rriT = rn'rj^, then sign (T) = sign (T'). 

To obtain the desired cancellation-free formula we therefore only need to compute the 
number of trek systems T' with the same trek system monomial as T. For this we 
introduce combinatorial objects that we call up-down cycles in T. Rather than giving a 
formal definition we describe their construction; see Figure |2] for an example. Let xq — > y 
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Figure 2. Solid directed edges denote parts of left paths in T, and dashed 
directed edges denote parts of right paths in T. On the left an up-down 
cycle is created, on the right, we traverse a right path up without hitting a 
further left path. 

be an edge in some left path Pl of T, and assume that xq — )■ y is not contained in any 
right path of T. Follow the directed path Pl downwards, i.e. with the direction of the 
graph, starting with xq — > y; since T has no sided intersection, we will never intersect 
another left path. If no vertex on P^ after the edge Xq ^ y is in any right path, then 
Xq ^ y is not in any up-down cycle. Otherwise, let yo be the first vertex on P]^ after 
y (possibly yQ equals y) that is also on some (necessarily unique) right path Pjj of a 
trek {P'^,P'pi) G T. Then follow P]j upwards, i.e. against the direction of the graph, from 
yo- If we never intersect another vertex on any left path (different from yo), then again 
xo — 2/ is not in any up-down cycle. Otherwise, let xi be the first vertex encountered 
when traversing upwards that is contained in any left path from T; say P'^. Then 
follow P'l down from xi, etc. Continuing in this manner we construct a sequence of non- 
empty paths Xo yo Xi ^ . . . in G that are alternatingly contained in left paths and 
right paths of T and that do not contain edges shared by some left path with some right 
path (because we always turn at the first possible vertex in each step). Only three things 
can happen: we traverse a left path downward to its final vertex without hitting a further 
right path, we traverse a right path upward to its initial vertex without hitting a further 
left path, or we eventually follow the original path Pl from somewhere above the edge 
xo — )■ y. In the final case, the sequence of paths closes up and the union of their edge 
sets is what we call an up-down cycle in T. The set of all up-down cycles in is denoted 
UD(T). For a more formal definition see the next subsection on the general directed case. 

Theorem 3.4 (Power-of-two phenomenon for acyclic digraphs). Let G = {V,D) be an 
acyclic digraph and let A and B he subsets of V of the same cardinality, and assume we 



have fixed linear orderings of A and B, as in Theorem 3j3_. Let T he a trek system in G 
without sided intersection. Then the number of trek systems T' G T{A, B) with mx' = rriT 
eg^/afe2lUD(T)l. 

A consequence of the previous two theorems is the following formula. 

Corollary 3.5 (Cancellation-free formula for acyclic digraphs). Let G = {V,D), A, and 
B be as in Theorem \3.3[ Then 

det Ea,b = ^^S^ (T)2l™mimT 

[T]^er{A,B)/~ 
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where the sum runs over equivalence classes of the relation ~ defined 6t/ T ~ T' if and 
only if mx = 'ht.t' ■ 

3.2. General directed graphs. For a general digraph G = {V, D) that may contain 



cycles, the analogues of Theorems |3.3| and |3.4| and Corollary |3.5| require some notions 
from [10]. 

Definition 3.6. Let A = {ai, . . . , a^} and B = {bi, . . . , b^} be /c-subsets of V with fixed 
linear orderings. A self-avoiding flow from A to i? in G is a pair F = (P, C) where P is 
a set of k self-avoiding and pairwise vertex-disjoint paths whose initial vertices exhaust 
A and whose final vertices exhaust B, and where C is a set of non-empty, self- avoiding 
directed cycles that are pairwise vertex-disjoint and that are also vertex-disjoint from all 
paths in P. The path component P gives rise to a bijection A ^ B, which through 
the linear orderings corresponds to an element tt G &k- The sign signF is defined as 
sign (tt) ■ (—1)1*^1, where |C| denotes the number of cycles in the collection C. 

In this definition a self-avoiding cycle is identified with its set of edges (without distin- 
guished initial vertex). We introduce a trek analogue of this notion as follows. 

Definition 3.7. Let A = {oi, . . . , a^} and B = {bi, . . . , bk} be /c-subsets of V with fixed 
linear orderings. A self-avoiding trek flow in the digraph G is a pair T = {F^, Ffj) where 
Fl is a self-avoiding fiow from a (necessarily unique) /c-subset S' of V to A, and F/j is a 
self- avoiding fiow from the same set S to B. We call S the set of tops of the trek fiow. 
We write T{A, B) for the set of self-avoiding trek fiows from A to B. 

Note that our notation T and T{A, B) is consistent with the acyclic case: there the 
cycle components Cl and Cr of Fl and F/j are empty, and the condition that the path 
component of F/, (respectively, F^) consists of self-avoiding and pairwise vertex-disjoint 
paths is equivalent to the condition that a trek system has no sided intersection. 

Definition 3.8. The sign of a self-avoiding trek fiow T = (F^, F^) from AtoB (relative 
to fixed linear orderings on these sets) is defined as sign (T) = signF^j ■ signF/,. (This 
depends on the orderings of A and B, but not on the linear ordering of the set of tops, as 
long as the same ordering is used for determining the signs of both F^j and F^). 

The trek flow monomial rriT is defined as the product of the variables ua for i & S, the 
variables Xij corresponding to all edges i ^ j used by paths and cycles in F^ and F/j, 
taken with the appropriate multiplicities. 

Note that, in particular, the fact that T is self-avoiding implies that each Xij appears 
with degree at most 2 in niT- Using a theorem by the third author [in], we will write 
det Tja,b as a rational function where both the numerator and the denominator are signed 
sums over self-avoiding trek fiows; see below. The following theorem says that no cancel- 
lation occurs. 

Theorem 3.9 (Positivity for general digraphs). Let G = {V,D) be a digraph and let A 
and B be subsets of V of the same cardinality, both with fixed linear orderings. If two 
self-avoiding trek flows T and T' from A to B satisfy rriT = m'rp, then sign (T) = sign (T'). 
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Next we compute the number of trek flows T' with the same trek flow monomial as T. 
Again, our main combinatorial tools are up-down cycles, which we define more formally 
at this point. 

Definition 3.10. Let T = (Fl,Fji) be a self-avoiding trek flow, and let E be the set of 

edges used by F^, or by but not by both. Define a directed graph T with vertex set E 
as follows: an element e = x ^ y E E used hj Fl has at most one outgoing arrow in F. 
If y is visited by some path or cycle in and has an incoming edge f = z ^ y there, 
then the arrow from e points to /. If y is not visited by any path or cycle in F/j (not 
even by an empty path based at y) and has an outgoing edge f = y ^ z in Fl, then the 
arrow from e points to /. In all other cases, e has no outgoing arrow in F. Similarly, an 
element e = x y E E used by F^ has at most one outgoing arrow. If x is visited by 
Fl and has an outgoing edge f = x ^ z there, then the arrow from e points to /. If x is 
not visited by F^ and has an incoming arrow f = z ^ x in F^, then the arrow from e 
points to /. In all other cases, e has no outgoing arrow. 

Now an up-down cycle in T is the support F C E C D of any non-empty directed cycle 
in the digraph F. We write UD(T) for the set of all up-down cycles in T. 

In fact, the definition of F forces every element of E to have at most one incoming 
arrow, as well. As a consequence, F is a union of vertex-disjoint directed cycles and 
paths. This implies that distinct up-down cycles in T do not share edges, a fact that we 
will use later. It is easily verified that this definition of up-down cycle agrees with the 
construction given in the acyclic case. One new aspect is that the edges of a directed 
cycle in the cycle component of F^ which avoid F^ (or vice versa) also form an up-down 
cycle. 

Theorem 3.11 (Power-of-two phenomenon for general digraphs). Let G = {V,D) be a 
digraph and let A and B he subsets of V of the same cardinality, both with fixed linear 



orderings, as in Theorem 3j_9. Let T be a self-avoiding trek flow in G. Then the number 



of trek flows T' G T(v4, B) with rriT^i = mx equals 21^°*^'^^ 

Finally, we state our general cancellation-free rational expression for det T,a,b- 

Corollary 3.12 (Cancellation-free formula for general digraphs). Let G = (V, D), A, and 
B be as in Theorem \3.!A Then 

. _ _ E[T].er(AW^^ign(T)2l^P(^)lmx 
^'"^ E[T].er(0,0)/^sign(T)2lUBmimT' 

where both run over equivalence classes of the relation ~ defined 6y T ~ T' if and only if 

Note that every self-avoiding trek flow from the empty set to itself is of the form 
((0, Cl), (0, Cr))- In particular, when G is acychc and there are no cycles, the denomina- 
tor contains only one term, corresponding to = Cr = 0, and our formula specializes 



to the formula in Corollary 3.5 
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4. Determinants of path matrices for cyclic graphs 



In this section we recall a generalization of the Lindstrom-Gessel-Viennot Lemma [7j 
to arbitrary graphs; for a different generalization see j2j. Thus let G = {V,D) be an 
arbitrary finite directed graph, and associate an edge variable Xij to each i ^ j E D. Let 
A = {ai, . . . , ak} and B = {bi, . . . , bk} be fc-subsets of V, and let M = Ma,b denote their 
weighted path matrix, i.e., the k x fc-matrix over M[[Aij | i — )■ j G D\\ whose (i, j)-entry 
is the sum over all directed paths in G from a, to bj of the product of the edge variables 
along the path. 

To express det(M) as a rational function in the edge variables, we recall from Section [s] 
the notion of self-avoiding flow from A to 5 in G, and we write J^a,b{G) for the (necessarily 
finite) set of all such flows. 

The determinant of the weighted path matrix from A to B is 
^fgTaMG) sign(F)mF 



Theorem 4.1 ([HI HO 

given by 



det(M 



A,B) 



EceJ-„(G) sign(C)mc' 
where the denominator is the sum over all self-avoiding flows consisting of vertex- disjoint 
cycles only. 

5. Proofs of the main results 



In this section we prove Theorems |3.9| and |3.11| and Corollary |3.12[ which immediately 
imply their acyclic special cases. Let G = (V, D) be a directed graph and let A = 
{ai, . . . , afc} and B = {bi, . . . , b^} be fc-subsets of V. We first deduce the corollary from 
the two theorems. 



Proof of Corollary \3. given Theorems \3.9\ and \3.11] We will apply Theorem |4. 1| with G 
replaced by the digraph H constructed from G as follows: in the disjoint union of the 
opposite graph with G add an arrow from any vertex of to the corresponding 
vertex i in G. These arrows get labels Uu, while every edge — )■ 2°^^ in G°pp gets the 
same label Aj,- as its opposite z — )■ j in G. With this notation, the submatrix E^^ b of our 



graphical model equals the path matrix Ma°pp b in Theorem 4.1 Similarly, self-avoiding 



trek flows in G from A to B are in bijection with self-avoiding flows in H from A"^^^ to 



B, and the bijection preserves signs. Thus Theorem yields the expression 

T^TeTiAB) sign (T)mT 



detS 



A,B 



ETer(0,0) sign (T)mT 



Theorem 3.9 shows that self-avoiding trek flows from A to B with the same trek flow 



monomial have the same sign, and Theorem 3.11 counts the number of such trek flows in 



each equivalence class. This gives the desired formula 

Em..rM,m/^sign(T)2lui^m 



detS 



A,B 



.er(0,0)/'- 



sign(T)2lUD(T)lr7iT 



□ 
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The remainder of this section focuses on the proofs of Theorems 3.9 and 3.11 The 



proofs are closely intertwined and mostly contained in the main body of the text. 

First we want to reduce our arguments to the case where A and B are disjoint. For 
this, we first modify G as follows. For each Oj, 6j, introduce copies a'-, b[ with edges aj — )■ 
a^, bi — )■ b[ and call the resulting graph G'. Let A' = {a'l, . . . , a'j,} and B' = {b[, . . . , 6'^} 
inherit their linear orderings from A and B respectively. Adding the new edges gives a 
bijection from self-avoiding trek flows in G from A to S to self-avoiding trek flows in G' 
from A' to B'. For Theorem 3^ we observe that this bijection preserves signs, and for 
Theorem 3.11 we observe that it also preserves the number of up-down cycles — indeed, 
the new edges are not part of any up-down cycle. Thus, we may drop the accents and 
write G,A,B,ai,bi for G',A',B',a'-,b'-, or in other words, assume A and B are disjoint 
subsets of V. 

Next fix a self-avoiding trek flow T = (F^ = (P^, Cl), Fr = (Pr, Cr)) from A to B 
with set of tops S = {si, . . . , Sk}, and choose the linear orderings such that Pl and Pr 
connect Si, i E [k] to and to bi, respectively, so that the path components of both 
flows have sign 1. We want to show that any self- avoiding trek flow T' from A to B with 
the same trek flow monomial as T has the same sign as T, and that the number of such 
self-avoiding trek flows is 2^^^^'^^^. 

To this end, color the edges of G with if they are in some path or cycle of F^, and with 
1 if they are used by F^. Delete edges of G without any color; these are irrelevant for what 
follows. Then contract all bi-colored edges in G. This may yield some additional isolated 
vertices (corresponding to fully contracted cycles all of whose edges are bi-colored); delete 
all isolated vertices. Contracting bi-colored edges and deleting resulting empty cycles and 
other isolated vertices gives a bijection between self-avoiding trek flows T' from A to B 
in the original graph satisfying mx' = and self-avoiding trek flows in the resulting 
graph having full support (i.e., containing all edges exactly once) and having S as set of 
tops. This bijection is sign-preserving because it does not affect the bijections A B 
and because it deletes cycles in pairs. Also note that the number of up-down cycles in 
T does not change, since up-down cycles by definition do not contain bi-colored edges. 
Again, we change notation so that G stands for the new graph. 

Next, consider any monochromatic cycle in G that has no vertices in common with 
paths or cycles of the other color. This cycle can be moved freely between the left and 
right cycle components of a self-avoiding trek flow with full support without changing 
the sign or the trek flow monomial. Moving such a cycle gives a factor 2 when counting 
self-avoiding trek flows of full support with 5* as set of tops, and this factor is accounted 
for by the fact that the cycle counts as an up-down cycle. So it suffices to prove our 
theorems for the simplified graph where all such cycles (edges and vertices) are deleted 
from G. Update G again. 

Finally, contract all vertices with exactly one incoming edge and one outgoing edge 
(necessarily of the same color). This can now be done without contracting a monochro- 
matic cycle (since such cycles always intersect a cycle of the other color after we have 
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removed the monochromatic cycles from the preceding according to the preceding para- 
graph). Again, there is a sign-preserving bijection, and the number of up-down cycles is 
invariant. 

After these modifications, the vertices in S have exactly two outgoing edges, one colored 
and one colored 1, and no incoming edges, the vertices in AUB have exactly one incoming 
edge colored (for A) or 1 (for B) and no outgoing edges, and all other vertices in G have 
exactly two incoming edges (of different colors) and exactly two outgoing edges (again of 
different colors). The resulting graph, again denoted G = {V, D), does not have any loops 
(but of course it may have directed cycles). 

We can recover T from the coloring by taking for F/j the 0-colored edges and for F^, the 
1-colored edges. We may therefore identify T with the Z/2Z-coloring. More generally, self- 
avoiding trek flows T' from A to B with the same monomial as T correspond bijectively 
to edge colorings D — t- Z/2Z with the same properties as T, namely: the incoming edges 
into A are colored 0, the incoming edges into B are colored 1, the outgoing edges of each 
vertex mV—{AUB) are colored with distinct colors, and so are the two incoming edges 
of each vertex in V \ {S U AU B). We identify such T' with their colorings. Then the 
sum T' + T (modulo 2) is an edge coloring where the sum of the colors of the edges 
entering any vertex is zero, and so is the sum of the colors leaving any vertex. Conversely, 
adding a coloring h with these properties to T yields a self-avoiding trek flow T' with the 
same monomial as T. Since the conditions on h define a subgroup K{T) of (Z/2Z)-^ the 
number of self- avoiding trek flows T' with the same monomial as T is a power of 2. The 



following proposition completes the proof of Theorem 3.11 

Proposition 5.1. To each up- down cycle C in T assign an edge coloring D — t- Z/2Z 
that colors the edges in C with 1 and the remaining edges with 0. The resulting up-down- 
colorings form a basis of K{T) as a vector space over Z/2Z. 

Proof. Let h he a. non-zero element of K{T) and let Xq — > yo be an edge contained in the 
left part of T that is colored 1 by h. Then yQ has a unique other incoming edge Xi — )■ yo 
that must also be colored 1 by h, and this edge must be in the right part of T. Then xi 
has a unique other outgoing edge xi — )■ y2 that must also be colored 1 by h, and this edge 
must be in T^, etc. All these edges are distinct, until you get back to xq — )■ yi. This gives 
an up-down cycle whose coloring can be subtracted from h to get a coloring with fewer 
edges colored 1. This proves that up-down colorings span K{T). Since distinct up-down 
cycles are edge-disjoint, up-down colorings are also linearly independent. □ 

To prove that the signs of T' and T are the same, it now suffices to prove that adding to 
T a coloring supported on an up-down cycle does not change the sign. For this we prove 



Lemma 5.3 To do this we first need to give a definition of signs of other combinatorial 



objects, which are used in the proof. 

Definition 5.2. Let 6 be an element of &k, let X,Y be subsets of [k] of the same 
cardinality, and let t; be a bijection from [k] — Y to [k] — X. We define the sign of the 
pair {S,v) as follows. Construct a directed bipartite graph H on two copies [k]i, [k]2 of 
[k] with 6 prescribing the arrows from [k]i to [A;] 2 ("down") and v prescribing the arrows 
from [k]2 ("up"). The graph H defines a bijection vr : X — )■ F which maps x G X to the 
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endpoint of the path in H starting at X. The sign sign (5, v) is defined as the sign of vr 
(relative to the natural linear orderings on X, F C \k]) times (—1) raised to the number 
of directed cycles in H. 

Lemma 5.3. Let v be as in the definition. Then sign {6, i')sign (5) depends only on v and 
not on 5 & &k- 

Proof. We prove the lemma by induction on A; — |F|(= k — \X\). It is trivially true if X and 
Y are both equal to [k]: then the graph H does not have any upward arrows and ir — 5, 
and sign {S, v) — sign (tt) = sign (S), so that sign {S, v)/sign (S) equals 1, independently of 
S. 

Now suppose that the lemma is true for all v : [k] — Y ^ [k] — X. Pick such an v, let 
X E X and y E Y, and let v' be the extension of v to [k] — Y + y sending y to x. Let 
5 E &k and construct H and tt as in the definition (from the pair 5,v). To analyze how 
sign {S,v') differs from sign {S,v) we distinguish two cases: 

(1) 7r{x) — y: In this case, adding the upward arrow y ^ x to H creates an additional 
cycle in H, and this contributes a factor —1 to the sign. Furthermore, the new 
bijection tt' : X — x ^ Y — y is the restriction of n to X — x. Hence we have 

sign (7r')/sign (tt) = sign ( JJ {i - x) {%{{)- yU = {-l)\ii^^\'<^}\+\{jeY\j<y}\_ 

\iex-x / 

In conclusion, sign {5, v')/sign {5, v) is this latter expression times —1 (for the ad- 
ditional cycle created in H). 

(2) tt{x) = z ^ y: In this case, no new cycle is created when adding to H the upward 
arrow y ^ x. Let u := n^^{y), so that in the new graph the path starting 
at u E X passes through u, y, x, z in that order. Then the new permutation 
tt' : X — X ^ Y — y equals tt on X — u — x and sends u to z. Hence we have 

sign(7r') . / n.eA-.-.(' -'0(^(0- --) 

sign 



sign(7r) \{u - x){y - z)Yl-^x_^_^[{i - u){7r{i) - y){i - x){7r{i) - z)] 

1 



sign 



= —sign 



{u -x){y- z) Yii^x-u-xlW) -y)ii- ^)] 

1 



^ _(_iy{«ex|i<a;}|+|{iey|i<jy}|. 

We conclude that the change of sign is the same in both cases, and that it does not depend 
on 5. □ 

We will use the following direct consequence of the lemma. 

Lemma 5.4. Fix a natural number k. Let So,Si be elements of &k and let Vo,Vi be 
bijections between subsets [k] — Yq, [k] — Yi and [k] — Xq, [k] — Xi, respectively. Hence 
— 1^1 for i — 0,1; we do not require that \Xo\ — \Xi\. Then we have the following 
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identity: 

sign (5o, vo) ■ sign {5i, Vi) = sign i;o) ■ sign {Sq, Vi). 
Now we return to our proof of positivity; the following arguments are clarified in Ex- 



ample 5^ below. Suppose that T' and T differ by an up-down cycle h, considered as an 
edge coloring as above. Label the vertices with incoming edges of /i-color and outgoing 
edges of /i-color 1 by the numbers 1, . . . ,k (we identify these with a first copy [k]i of [k]), 
and similarly label the vertices with outgoing edges of /i-color and incoming /i-color 1 by 
1, . . . ,k (identifying them with a second copy [k]2)- Note that there are, indeed, equally 
many of both, since the edges of /i-color 1 and T-color form k vertex-disjoint paths from 
the first set to the second set. Denote the corresponding bijection hj Sq ^ &k- Similarly, 
the edges of /i-color 1 and T-color 1 give a bijection Si E &k- Now some directed paths 
with T-color emanating from [k]2 (hence with /i-color 0) may "loop back" to [A;]i; this 
gives an bijection vq from a subset [k]2 — Yq of [k]2 to a — Xq. Similarly, we obtain an 
injective map vi from a subset of [A;] 2 — Yi into [k]i — Xi by following directed paths of 
T-color 1 emanating from [A;] 2. 

To change T into T' the colors in the up-down cycle are interchanged. This means 
that the roles of 60 and 61 are interchanged. The total change in trek flow sign is exactly 
the change from sign (5o, ^o)sign (5i, into sign (5i, i;o)sign (5o, t'l)- Indeed, the change 
in sign of the left flow is exactly sign (5i, i;o)/sign (5o, i^o), and the change in sign of 



the right flow F/j is exactly sign (5o, ^1) /sign ((5i, Now Lemma 5.4 implies that that 
T and T' have the same sign. This concludes our proof of Positivity for general digraphs. 
We conclude the paper with an example illustrating the arguments just given. 



Example 5.5. In Figure [3]we have k = 3, the elements of [3]i are denoted 1, 2, 3, and the 
elements of [3]2 are denoted l',2',3'. On the left is a fragment of T with 0-colored edges 
solid and 1-colored edges dashed. On the right is the corresponding fragment of T'. 

The bijection 60 maps 1 h-> 1', 2 2', 3 (-)■ 3' and is depicted by straight solid arrows 
on the left and by dashed solid arrows on the right. The bijection 61 maps 1 1— )■ 2', 2 1— )■ 
l',3 I—)- 3' and is depicted by straight dashed arrows on the left and by straight solid 
arrows on the right. Together these form the up-down cycle. The partial map up vo maps 
2' H-)> 3 (curved solid arrow) and the partial map vi maps 1' 1— )■ 2, 3' 3 (curved dashed 
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arrows). The remaining arrows are fragments of paths leading from the set of tops to the 
up-down cycle and from the up-down cycle to AU B. 



The map ttq constructed from (So^Vq) as in Definition 5.2 can be read off from the 
left-hand picture: it maps 1 t— t- 1', 2 i— )■ 3' and has sign 1; the number of cycles in the 
digraph Hq as in that definition is zero. So sign ((5o, fo) = 1- Similarly, tti constructed 
from vi) is 1 I— !■ 2', with sign 1, and Hi has two cycles 21' and 33', so sign {6i,vi) = 1, 
as well. 

On the right, the roles of 6o and 6i are reversed. The pair {6i,vo) leads to no cycles 
and to the map 1 3', 2 t-i- 1' with sign —1. Thus sign t;o) = —1. Similarly, {6o,vi) 
leads to the map 1 i— )• 2' with sign 1 and to one cycle 33', so sign(5o,'yi) = —1- Note 



that sign (5o, 'Uo)sign Vi) = sign {6i, Vo)sign (5o, Vi) in accordance with Lemma 5.4, and 
this shows that the sign of T does not change under swapping the and 1-labels in the 
up-down cycle, resulting in T'. 
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